Quantum chaos in an ion trap: the delta-kicked harmonic oscillator 
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We propose an experimental configuration, within an ion 
trap, by which a quantum mechanical delta-kicked harmonic 
oscillator could be realized, and investigated. We show how 
to directly measure the sensitivity of the ion motion to small 
variations in the external parameters. 

PACS: 05.45.+b, 03.65.Bz, 42.50.Vk 



In classical mechanics, deterministic chaos is often 
most simply described as exponential sensitivity to initial 
conditions, meaning that initially neighbouring classical 
trajectories diverge extremely rapidly with time. Due to 
the necessity of preserving the inner product, this kind 
of divergence between two possible initial states cannot 
occur quantum mechanically. The question of what then 
constitutes the quantum mechanical equivalent of chaos 
immediately arises. An interesting proposal by Peres |l]] 
is to examine the initial state \tp} evolving under two 
slightly differing (classically chaotic) Hamiltonians Hi 
and H 2 . Defining Ui_ 2 (t) as the corresponding unitary 
evolution operators, the overlap 
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is predicted to behave very differently, depending on 
whether the initial state is in a stable or chaotic area of 
phase space. Thus O is a measure to distinguish between 
regular and irregular quantum dynamics. In fact much 
work on the subject of quantum chaos has been carried 
out theoretically; experimental realizations however re- 
main somewhat scarce although there have recently 
been pioneering successes in atom optics || and in meso- 
scopic solid state systems j|. In this Letter we propose 
a realizable experimental configuration, with which one 
can measure O directly. The system proposed is a single 
ion trapped in a harmonic potential, subject to periodic 
kicks from a standing wave laser. This is a quantum 
delta-kicked harmonic oscillator, a system capable classi- 
cally of stochastic dynamics, including Arnol'd diffusion 
H under certain resonance conditions || . 

Trapped ions are in many ways an ideal choice of sys- 
tem for the study of fundamental aspects of quantum me- 
chanics. One can take advantage of the small dissipation 
in this system, together with the possibility of coherent 
manipulation of the ion's motional state. Trapped ions 
have been used in recent experimental demonstrations 
of the generation of non-classical states of motion 0, 
quantum logic gates and tomography of the density 
matrix |)J. In addition there have been theoretical pro- 
posals for investigation of localization J10| and irregular 



collapse and revival dynamics |TT| due to quantum chaos 
in this system. 

In this Letter we will first describe a general proce- 
dure for determining O as a function of time. We then 
describe explicitly a particular system, the delta-kicked 
harmonic oscillator, how it may be implemented within 
an ion trap, and how to carry out our general procedure 
for determining O. Finally we display some numerical 
results, showing what one would expect to see when car- 
rying out such an experiment. Markedly different results 
are indeed observed numerically, dependent on whether 
the initial condition is in a classically stable or chaotic 
area of phase space. 

We first consider a general Hamiltonian of the form 
H = Hi\gi)(gi\ + H 2 \g 2 )(g 2 \, where \gi) and \g 2 ) are 
stable electronic ground states of a single trapped ion. 
The state of the ion is set initially to be: 
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Here \a) and \j3) are states of the ion's motion. In the 
applications described in this Letter, \a) and \(3) are 
coherent states; IV'(O)) is then a Schrodinger cat state, 
which has been achieved experimentally for a trapped 
ion 0. After a time t, the initial state evolves to 
\m\ = [\9i)Ui(t)\<x) + |<?2>?7 2 (t)|/3>]/V2, where U x {t) 
and U 2 (t) are the time evolution operators derived from 
Hi and H 2 , respectively. A 7r/2 pulse is applied (Ramsey 
type experiment) to the ion, yielding 



\m)' = \{\9i) [Ui(t)\a)-U 2 (t)\0) 
\g 2 ) [Ui(t)\a) + U 2 (t)\f3)]}. 
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The probability for the ion to be in state \gi) is thus 
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Similarly, if we set |V>(0)) = (\gi)\a) + i\g 2 )\/3))/V2, the 
corresponding final probability is given by 
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\Ul(t)Ui(t)\a)]} . (5) 



By determining P g and P' g , one can clearly deduce 

\(P\Ul(t)Ui(t)\a)\ 2 . If Hi and H 2 are slightly differing 
chaotic Hamiltonians, and \a) = |/3), then we have O, as 
defined in Eq. ([[]). We also note that where H 2 is the 
simple harmonic oscillator Hamiltonian, this collapses to 
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2ttQ((3), where Q(f3) is the Q function for various initial 
(3 of the pure state Ui(t)\a). By repeated measurements 
one can therefore determine the Q function's evolution 
in time |l2j . 

Our proposed model system is a harmonic oscillator 



Hn = 
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periodically perturbed by nonlinearly position-dependent 
delta-kicks; 



Hi = K cos(fcx) — nr ) 



(7) 



so that H = Hq + Hi. Here x is the position, p 
the momentum, m the mass, v the oscillator frequency 
k = 2tt/X the wavenumber, t — time, r the time delay 
between the kicks, and K the kick strength. Under the 
resonance condition vt = 2Trr/q (r/q is a positive ratio- 
nal, where q > 2), classically there are thin channels of 
chaotic dynamics in the phase space The resulting 
Arnol'd stochastic web [see Fig. 0(a]j spreads through 
all of phase space; Arnol'd diffusion || can occur in sys- 
tems of less than two dimensions when the conditions for 
the KAM (Kolmogorov, Arnol'd, Moser) theorem [jl3| are 
not fulfilled, as is the case here 0. The corresponding 
quantum mechanical system has also been studied theo- 
retically Q [see Figs. |l](b,c,d) for the time averaged Q 
function of this system]. 

To construct such a system quantum mechanically, 
which can also be used to carry out the procedure de- 
scribed in Eqs. we begin with a single ion in a 
harmonic potential (e.g. a linear ion trap |7|]); in addi- 
tion we require a time dependent standing wave laser 
configuration. The ion has two ground states and two 
excited states, and the laser is elliptically polarized [see 
Fig. ||(a)]. The <r+ and cr_ polarized contributions thus 
separately couple two different two level systems, with 
different Rabi frequencies: 
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where ojq is the transition frequency between the elec- 
tronic states \ej) and \gj), u>l is the laser frequency, and 
^1,2 (t) are the (time dependent) Rabi frequencies. In a 
rotating frame defined by U = exp[— iujjj, Y^j=i(\ e j) ( e j l — 
\gj)(gj\)/2], and in the limit of large detuning |A| = 
\u>l — (jJq\ ^> \Qi t 2(t)\, |ei) and |e2) can be adiabatically 
eliminated to give: 

h 2 

^ = ^° + 8A ^ Qj ^ 2 [C ° S {2ki) + 1] l9j){9jl (9) 
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The laser is rapidly and periodically switched, giving a 
series of short Gaussian pulses: 



(10) 



which approximate a series of delta kicks in the limit a — > 
0. Note also that we require a ^ 1/A, otherwise the laser 
is too spectrally broad, making adiabatic elimination of 
|ei) and |e2) impossible. Thus, finally we have: 

2 oo 

H = H + Kj [cos (2kx) + l]\ 9j )( 9j \ S(t-nr), 
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(11) 

which corresponds almost exactly to Eqs. (fU?]), for two 
different Kj = hay/nQ 2 /8A. There are extra \gj)(gj\ 
terms, but these will only contribute phases to the evo- 
lution of the initial state of Eq. (||), and can easily be 
accounted for. 

Taking |a) = \/3), the initial state of Eq. (|^) thus 
evolves as 



|^(nr)) = l^e-- ^ F?\ 9j )\a) 



(12) 



where the Floquet time evolution operators Fj are given 
by: 



Fa 



-iaJ ai/T —iKj cos[2?7(d^+d)] \gj) {gj\/ v^r/ 2 



(13) 



where a) and a respectively create and annihilate a single 
phonon quantum, and the common phase term e r lVT I 2 
has been dropped. The Kj = r^?7 2 cr-\/27r/8A are dimen- 
sionless kick parameters, which, with vt, determine fully 
the phase space behaviour of the classical delta-kicked 
harmonic oscillator || . In the quantum mechanical prob- 
lem there is an additional parameter, the Lamb-Dicke 
parameter r) = ky / h/2mv. As rj 2 oc Ti, by progressively 
reducing r\, one can explore the transition from quan- 
tum to classical chaos |y| . This can be accomplished by 
"tightening" or "loosening" the trapping potential, i.e. 
increasing or decreasing the trapping frequency v. 

After n kicks, we perform a 7r/2 pulse between the lev- 
els |<7i) and \g"ih, by e.g. a Raman transition or a magnetic 
field [see Fig. g(b)]. By fluorescence, using an auxiliary 
level |/), with repeated measurements one can determine 
P g and P' g [see Fig. |(c)], as defined in Eqs. (§||). 
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where Sk = Ki — K\. From Eq. ([TJJ) one can easily extract 
the overlap O = \{a\F^ n F^\a)\ 2 . 
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In order to relate the quantum mechanical behavior of 
the system to the classical one, we have to use an initial 
condition equivalent to the classical x(0) and p(0). We 
use a coherent state \a), where a can be expressed as a — 
[kx(0) + ikp(0) / mv] / 4r] . Thus we can see that when rj is 
small, a is large, and \a) is therefore more macroscopic, in 
some sense more classical. This can be seen by comparing 
Fig. |l](c) with Fig. [j](d); for r/ = 0.5 population "tunnels" 
through a classically forbidden area, which does not occur 
when ry = 0.25 @. 

Figure || shows the values of P g and P' g that one would 
measure for this scheme, and the value of O that one 
would thus obtain, for r\ = 0.5 after 0-1000 kicks. The 
plots obtained are clearly different, depending on whether 
the initial condition is classically unstable, as in Figs. 
^(a,b) where the corresponding classical initial condition 
is a hyperbolic fixed point, or stable [Figs. ||(c,d), elliptic 
fixed point]. This is already noticeable in the plots of 
P g and P' g) before O is extracted [Figs. ||(a,c)]. In line 
with previous numerical work for the kicked top pj], O 
decays for an unstable initial condition, and undergoes 
quasistable oscillations for a stable initial condition. 

As O is a measure of how close the two parallel evo- 
lutions are at a given time, it can be seen that if the 
initial condition is classically unstable [Fig. ||(b)] the two 
states become rapidly increasingly orthogonal (more "far 
apart"), whereas in the case of a stable initial condition, 
for some time the difference between the states remains 
on average about the same. This in some sense corre- 
sponds to the classical definition of chaos, where under 
the influence of the same dynamics, very slightly differ- 
ent initial states diverge rapidly if their origin is in an 
unstable area of phase space jjj . 

Figure || shows the same for n = 0.25. In line with the 
fact that this is more in the semiclassical regime than 
Fig. |[ the decay [Figs. f|(a,b)] is more rapid, and the 
oscillations [Figs. §|c,d)] are more stable. The slow decay 
of the quasistable oscillations when r) = 0.5 [Fig. ||(d)] 
can be traced back to the tunneling that takes place in 
this regime (see Fig. |l|), absent when i] = 0.25. 

Numerically the procedure is carried out in a trun- 
cated Fock basis of 400 states when 77 = 0.5, or 800 when 
1] = 0.25. Increasing the size of the Fock basis does not 
qualitatively change the observed dynamics. 

In conclusion we have shown a general procedure for 
determining the overlap parameter O originally proposed 
by Peres Q. We have described explicitly how O could 
be determined for the delta-kicked harmonic oscillator, a 
classically chaotic system. We have described how a sin- 
gle ion trapped in a harmonic potential could be a practi- 
cal experimental realization of the delta-kicked harmonic 
oscillator, and how our scheme for determining O is re- 
alized in this configuration. In particular, our scheme 
presents a direct way for determining O, by virtue of the 
fact that we effectively have two Hamiltonians running 
in parallel, within the same experimental system. 
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standing light wave see ||. 

The tunneling effect can be measured by tomographic 
methods (see [|s|,|l2|). One can also also observe the tun- 
neling to the center of the web using the present scheme 
with (3 — in Eq. (^) and Ui(t) the free evolution of the 
harmonic oscillator [compare Eqs. (If) and (at)]. 
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FIG. 3. (a) P g (solid line) and P' g (dotted line) for 
the initial condition a — n/2r] (rj = 0.5, «i = 0.2, 
k 2 = 0.225). (b) O (solid line), Re((Q|F" t F; 1 |q)) (dashed 
line), and lm({a]F" t F"]a)) (dotted line) for the same ini- 
tial condition, (c) and d), same as a) and b), for the initial 
condition a — in/ri\f3. 



FIG. 1. Stroboscopic Phase space portraits (from kick 
to kick) for the delta-kicked harmonic oscillator, where 
r/q — 1/6 and k — 0.2, after averaging over 11000 kicks. 
The position x and momentum p axes are in units of A and 
mvX respectively, (a) Poincare surface of section, showing the 
classical stochastic web. Regular dynamics take place within 
the cells defined by the web. (b) Time averaged Q function 
for the (unstable) initial condition \a) where a = n/2rj and 
r\ = 0.25, centred at (1,0). (c) Time averaged Q function 
for the (stable) initial condition a — iiv/r]\/3 ((0, 2/\/3)) and 
r\ — 0.5. (d) As for (c), where r\ = 0.25. Note the "tunneling" 
out of the original ring of cells for the larger value of r\. 
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FIG. 4. Corresponds exactly to Fig. j^, except that 
r) = 0.25 
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FIG. 2. Atomic level configuration proposed for our proce- 
dure, (a) The atom first experiences a series of short laser 
pulses from an elliptically polarized standing wave. Levels 
|<7i) and |ei) are coupled by cr+polarized light, and levels \g2) 
and |e) by <r_ polarized light, where the corresponding inten- 
sities differ, (b) After a definite number of kicks, the levels 
\gi) and | experience a 7r/2 pulse (e.g. a Raman transition, 
using the auxiliary level \r)). (c) The population of |<ji) is 
determined by fluorescence, using the auxiliary level |/}. 
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